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ABSTRACT 



We present an early result from an automated search of Kepler eclipsing binary systems for circumbinary companions. An intriguing 
tertiary signal has been discovered in the short period eclipsing binary KIC002856960. This third body leads to transit-like features in 
the light curve occurring every 204.2 days, while the two other components of the system display eclipses on a 6.2 hour period. The 
variations due to the tertiary body last for a duration of -1.26 days, or 4.9 binary orbital periods. During each crossing of the binary 
orbit with the tertiary body, multiple individual transits are observed as the close binary stars repeatedly move in and out of alignment 
with the tertiary object. We are at this stage unable to distinguish between a planetary companion to a close eclipsing binary, or a 
hierarchical triply eclipsing system of three stars. Both possibilities are explored, and the light curves presented. 

Key words. stars:individual: KIC002856960, stars: general, binaries: close, binaries: eclipsing, planets and satellites: general 
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1. Introduction 

With approximately a third of main sequence stars in the galaxy 
found in binary systems (Lada 2006) but a much smaller per- 
centage of known planets in such, planets in binary systems rep- 
resent an important area of investigation. Transiting circumbi- 
nary planets can provide a wealth of information about them- 
selves and their host systems, from the extreme system flatness 
necessary for planetary transits and binary eclipses to be seen 
to potential absolute dynamical solutions once timing variations 
are taken into account. As such they are primary targets for at- 
tempted detection. The Kepler mission has already discovered 
three circumbinary planets ( |Doyle et al.|2011 Welsh et al.|2012| l 
and provides the opportunity for many more. The components of 
these three systems are all well separated, leaving the parameter 
space of transiting planets around contact and near contact bina- 
ries as yet unexplored. 

Triple stellar systems are also of great interest, and several 
have recently been detected through the Kepler data (Derekas 
|et al.|20lT Steffen et al.|201 1 1. If such systems transit, their sig- 
nals can be similar to those from planetary systems, and so a 
search for planets in binaries may also lead to detection of triple 
stars. These can provide insights into the stability, dynamics and 
evolution of complex systems. Recently such a triply eclips- 
ing triple was discovered and characterised using data from the 



Kepler mission ( Carter et al. 201 1| ), providing rare observational 
constraints for low mass stars with fully convective interiors. 

With the aim of expanding the parameter space of known cir- 
cumbinary planets a search was performed on the eclipsing bi- 
naries listed by the Kepler Eclipsing Binary Catalogue (KEBC) 
( |Prsa et aL 201 l|[Slawson et al.|201 l| l, initially on those systems 
not classified as detached. Here we present one early discov- 
ery from the search: KIC002856960 (RA(J2000) 19:29:31.52, 
DEC +38:04:35.9, Kepler magnitude 15.615), a system com- 



posed of either a close eclipsing binary orbited by a transiting 
planetary companion or low mass dwarf star, or a triply eclips- 
ing group of at least three low mass stars. This result and others 
in the future should prove of value to models of complex sys- 
tems, where gravitational interactions of several stars can lead 
to unusual perturbations (e.g. Quintana & Lissauer (2006) in the 
planetary case), and as such provide fundamental observational 
information to general formation models. 



2. Observations 

2.1. Data 

Observations of the KEBC systems were taken by the NASA 
Kepler satellite, a mission producing extremely high precision, 
near continuous light curves of -155,000 stars on the level of 
20 ppm ( |Koch et al.|2"0T0}|Batalha et al.|2010||Haas et al.|2010) . 
This mission began science operations on 13 May 2009, and data 
up to the end of Quarter 6 (22 Sep 2010) are publicly available 
on the NASA Data Archiv^] Due to technical constraints, the 
satellite must reorient itself each quarter year, and so data is pro- 
vided as a separate file per object and per quarter, each quarter 
consisting of three months of observations. Supplementary in- 
formation on targets and their designations are available from 
the Kepler Input Catalogue (KIC) ( |Brown et al.|201 1| >. The data 
presented were taken in long cadence mode (29.4 minute expo- 
sures) over a period of 498 days, and represent 2265 1 good ex- 
posures when small gaps in the data are taken into account. See 
the Data Characteristics Handbook on the NASA Data Archiv^] 
for precise observation periods. In raw format, the data show a 
variety of noise signals, the details of which are explained in the 

1 http://archive.stsci.edu/kepler/ 

2 http://archive.stsci.edu/kepler/documents.html 
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Short-Period Phased Data 




Fig. 1. The binary variation (eclipses of A and B by each other), 
shown as phase folded data from all 6 quarters (scale shown is 
relative offset to the median of each quarter). The tertiary event 
has not been removed, and as such some points are omitted be- 
low the figure window. These are displayed more clearly in Fig. 
[2 Error bars are not plotted, as the spread of points gives a much 
clearer indication of the true error. 



Kepler Data Processing Handbook 2 . The removal of the major- 
ity of the instrumental parts of these was performed using co- 
variance basis vectors (CBVs), enacted using the PyKe python 
modul^] This method attempts to remove instrumental effects 
by using variations observed in several stars neighbouring the 
object of interest, and as such should robustly maintain astro- 
physical signals. For a full explanation of the method see the 
software descriptioij^] 



2.2. Eclipsing Binary Signal Removal 

Data were then processed to remove the binary light signature 
from each quarter separately (hereafter 'whitening'). First, data 
were phase folded using the period provided by the KEBC. The 
phased curve was then binned using 200 equal width bins, and 
the median of each bin determined. As any points exhibiting ter- 
tiary transit events will be distributed across the light curve when 
phase folded on the binary period, taking the median will ex- 
clude them from the whitening process. The median of each bin 
is then subtracted from each point in the bin, resulting in a series 
of data points distributed around a median adjusted flux of zero. 
This method has the advantage of requiring no knowledge of the 
signal being removed beyond its rate of variance (required to al- 
low the bin number to be selected appropriately), and as such 
can be applied to many variable objects. The whitened quarters 
were then concatenated and searched for tertiary events. An au- 
tomated search method (to be described fully in Armstrong et al. 
2012b,in prep) combined with by-eye inspection of the results 
led to the discovery of the object presented here. 



3 http://keplergo.arc.nasa.gov/PyKE.shtml 

4 http://keplergo.arc.nasa.gov/ContributedSoftwareKepcotrend.shtml 



3. Results 

3.1. Overview 

Periodic signals separate from the previously known short pe- 
riod variations shown in Fig. [T] have been detected from the ob- 
ject KIC002856960. There are three new events associated with 
the long period signal seen in the public data, in quarters 1, 4 and 
6, and these are separated and presented in Fig. [2] They display 
a unique appearance, indicative of the presence of at least three 
bodies in the system. 

The short period variations have not been amenable to mod- 
elling by the catalogue team. As such, while we proceed under 
the assumption that the source of the short period variation is an 
eclipsing binary system (of a perhaps complex nature), no binary 
parameters beyond the ephemeris of the short-period variation 
are provided. We have derived an updated binary orbital period 
for this signal using all of the publicly available quarters ( 1 -6) 
via x 2 minimisation of each bin in the binned data (see Sect. |2.2)> 
and present this along with the catalogue value for To in Table 111 
The original binary signal is approximately 1% in depth, notably 
shallower than the newly detected signal. 

From here we shall refer to the bodies in this system as 
Components AB (for the two members of the binary system), 
and either b or C depending on whether the third body is con- 
sidered as a planet or star respectively. We limit ourselves to 
considerations of three bodies only, for simplicity. 

Each primary new event (hereafter 'crossing', as each event 
consists of a crossing with the binary orbit by a third body) con- 
tains multiple transits of A or B before C, or b/C before A or B 
(the distinction depends on the system architecture, see Sect. |4]i. 
Ephemerides for the three observed crossings are given, along 
with durations and depths, calculated from the deepest point ob- 
served in each crossing. The implications of these results for our 
interpretation of the nature of the system are discussed in Sect. 
|4] The crossing duration is given for comparison with the binary 
period - the first crossing in the quarter 1 light curve is cut off 
by the end of the quarter before completion of the crossing, and 
hence a duration is not provided. 

If secondary eclipse events in the long period signal (of com- 
ponent b/C by AB, or AB by C) are present, they are constrained 
by the variability of the data to have an upper limit to their depth 
of -0.4%, and no evidence for any is observed in the phased light 
curve. 

Time based parameter errors are dominated by the long ca- 
dence, in the absence of a full model. Crossings are taken as be- 
ginning midway between the data point immediately before flux 
variations are seen and its successor (the choice of these points 
is clear on viewing the data), and the absolute errors set such that 
the beginning of ingress could happen anywhere between these 
points. End times are derived similarly. All parameter errors de- 
rive from formal propagation from these values or the raw flux 
errors as appropriate. 



3.2. Blending 

It has been noted that the depths for each transit are markedly 
different; this is expected but quarter to quarter variations in 
the depth, along with the large size of the Kepler pixels of ~4 
arc sec (Koch et al. 2010| l can cause suspicion as to whether 
the object involved is blending with a background source. The 
photometric aperture used changes between quarters to opti- 
mise the photometric precision, meaning that a significant back- 
ground source may contribute a different amount of flux each 
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Fig. 2. The three tertiary crossings available in the public data, in quarters 1, 4 and 6. Crossings are referred to as Crossing 1, 2 and 
3, from left to right respectively. Light curves are displayed after removal of the binary signal. Dashed lines are simple connections 
of data points, for clarity, and do not stem from a model. 



Observed Aperture-Blending Analysis 




Table 1. Derived parameters for the KIC002856960 system 



Fig. 3. Images formed from the observed Q4 aperture for 



KIC002856960, using the method as utilised by |Torres et al. 
(2010 1. The direct image shows the location of the primary 
source, whereas the difference image shows the location of the 
flux variation. No discernible difference is observed. 



quarter. KIC002856960 has been subjected to a centroid analy- 
sis as utilised by |Torres et al. ( 2010| >. This involves creating a 
difference image across the pixel aperture between a selection 
of points in and out of transit. This difference image then shows 
the 'true' location of the flux variation. This is compared to a 
direct image formed from summing the same points. The direct 
image shows the location of the primary source of flux. Fig. [3] 
shows the difference image produced, and the comparison direct 
image. Blends originating from an offset of order 1 pixel are eas- 
ily ruled out by inspection. The fit of a pixel response function 
(a modified point spread function - Bryson et al. ( 2010| l) limited 
the region where a blending source could be located to within a 
radius of order 0.3 pixels from the target, but was complicated 
by the presence of the binary signal. Despite the difficulties this 
does rule out a large region of the photometric aperture as the 



Parameter 


Units 


Result 


Error (absolute) 


Period (AB) 


d 


0.2585082 




T (AB) 


BJD 


54964.652115 




Period (b/C) 


d 


204.2163 


0.030 


To (b/C) 


BJD 


54997.743313 


0.024667 


Duration (b/C) 


d 






Crossing 2 




1.2464 


0.0204 


Crossing 3 




1.2669 


0.0204 


Depth (b/C) 


% 






Crossing 1 




4.1171 


0.0588 


Crossing 2 




4.8179 


0.0669 


Crossing 3 




6.3207 


0.0749 



Notes. Parameters are labelled AB or b/C, standing for signals asso- 
ciated with the binary stars (AB) or companion (b/C). BJD stands for 
Barycentric Julian Date. T given as T - 2400000. A crossing 1 dura- 
tion is not provided as this crossing is cut off by the end of a quarter 
before completion. 



location of a background varying source. 

The nature of the tertiary signal requires that three bodies are 
present in the system; for it to arise from a separate system to the 
binary flux variations, while being so closely positioned, is ex- 
tremely unlikely. The depth of the tertiary signal also indicates 
that, if it was caused by a background source, the contaminating 
source would need to be of very significant brightness compared 
to the target. No evidence for such a source is seen. As such, we 
find it very likely that KIC002856960 is indeed the source of the 
signals presented in this paper. 



4. Discussion 

4.1. Summary 

Our primary interest lies in the tertiary signal detected. The pub- 
lic data contains three crossings (Fig. [2} appearing in quarters 1, 
4 and 6. These crossings occur on a period much longer than that 
of the binary AB, and are indicative of the presence of a com- 
panion third body b/C either crossing or being crossed by the 
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binary AB. The short period binary completes ~5 orbits during 
the course of the crossing, meaning that during a crossing multi- 
ple transits will be produced as the components A or B move in 
and out of alignment with the slower moving third body b/C. The 
shape of the transits seen within a crossing will be dependent on 
the velocity and relative positions of all components, as well as 
the components' shape and surface brightness. In the absence 
of a full model or complete mapping of the flux variations (due 
to the long cadence data only just resolving the inner transits 
within each crossing) we have presented the simplest observable 
parameters available: ephemeris, depth and overall crossing du- 
ration. Since the ingress/egress times are expected to depend on 
the location of all three bodies at the start/end of each tertiary 
event, an accurate period cannot be determined with only three 
crossing events. This situation will be largely improved once the 
currently private data becomes accessible. 

The structure and nature of this system remain in doubt. Here 
we discuss the two most likely system geometries: an object in 
orbit around a close binary pair (system (AB)b or (AB)C), or 
a close binary pair in orbit around a larger star (system C(AB), 
where C is the more massive star in this case). We consider these 
two options in turn, acknowledging that there may be other more 
unusual configurations possible. The equal depths of the AB bi- 
nary eclipses (Fig. [TJ imply that, if the binary components AB 
were on the main sequence, they should be of roughly equal 
mass. Given their short period, we find it plausible that the binary 
components may have evolved away from the main sequence, al- 
though without further data this cannot be constrained. In either 
case, the equal depth of the AB binary eclipses shows that the 
components A and B must have the same surface brightness. 



4.2. Circumbinary Object 

In this scenario the tertiary body b/C is of lower mass than the 
combined binary pair AB. A secondary eclipse of b/C by A or 
B that is not observable implies that the majority of the flux 
must be emitted from the components AB. Here we use the KIC 
colour derived parameters to explore the nature of this complex 
system. 

The KIC temperature (4733K) and stellar radius imply a star 
of luminosity O.26L . If both central stars AB were of this type, 
the 0.4% limit on secondary eclipses limits the luminosity of 
b/C to be ~0.002L o in this scenario. This is consistent with an 
object of spectral type M5 or cooler ( |Reid & Haw ley 2005 



Presenting a specific companion radius is beyond our capabi 
ity at this stage, but below we present example companion radii 
leading this configuration to support a circumbinary dwarf star, 
brown dwarf or perhaps transiting planet around a close binary 
pair, in the configuration (AB)b/C. The extremely shallow AB 
eclipses evident in the binary lightcurve could possibly be ex- 
plained as the result of a grazing geometry. 

The deepest observed transit (6.32%), on a single star of KIC 
stellar radius 0.757,R G ( |Brown e~a l. 2011), implies a compan- 
ion radius of 1.85/?/, or 0.19^ G . Assuming AB components of 
equal radii (and equal surface brightness as shown by the depth 
of their mutual eclipses), dilution would increase the companion 
radius to Q.21R Q . Given the very short period of AB it is possi- 
ble that the colour derived stellar radius will be an overestimate, 
allowing for the companion radius to decrease. There remains 
the possibility of the observed eclipses of A or B by component 
b/C being grazing, in which case the 'true' transit depth would 
be greater and the companion radius also. 



4.3. Circumstellar Binary 

This scenario assumes that the tertiary body is of greater mass 
than the binary pair. As such it would be the source of the ma- 
jority of the flux, and can be taken to be the object described by 
the KIC (consistent with a K star). This would then be orbited by 
the lower mass stellar binary pair AB, in the architecture C(AB). 
Taking the KIC stellar radius of the primary K star component 
C, we again derive a radius for a single component of the bi- 
nary system A or B of Q.19R Q . The system would then consist 
of a central K star orbited by a close pair of remarkably small 
dwarf stars. This has the advantage of explaining the shallow 
AB binary eclipse depths observed through dilution by the K 
star. The lack of secondary eclipses (eclipses of A or B by com- 
ponent C) imposes a problem however. To produce secondary 
eclipses of order 0.4% or less, such that they would not be vis- 
ible in the data as shown, requires that each component of the 
binary be contributing less than 0.4% of the system flux (for to- 
tal eclipses). This is not possible given the -1% binary eclipse 
depths observed (where the components would each need to be 
contributing at least 1% of the system flux). As such this con- 
figuration requires that the secondary eclipses be grazing, and 
thereby the primary eclipses also. Without further modelling we 
cannot take either configuration further. 

4.4. Conclusion 

At this stage we are unable to conclusively distinguish between 
the two system architectures (AB)b/C or C(AB) as outlined 
above. We also acknowledge that there may well be other sta- 
ble system configurations unexplored here. Spectral follow up 
and the currently private Kepler data, in combination with the 
presented lightcurves, will allow for a much greater understand- 
ing of this intriguing system. When its nature is fully under- 
stood the absolute determination of many of the system param- 
eters will be possible. This includes the companion mass in the 
(AB)b/C case, if transit timing variations as used in |Doyle et al.| 
(201 lj) and Welsh et al. (2012]) prove significant. We hope that 
this will prompt follow up studies, and will continue to develop 
our understanding of the object as more data becomes available. 
If the (AB)b case proves true, this system represents the first 
discovered transiting planet around such a relatively short pe- 
riod eclipsing binary. If the C(AB) or (AB)C cases prove true 
this system will contain stars with some of the lowest masses 
accurately determined. These will be low enough to imply fully 
convective interiors, an area where directly measured physical 
properties are few (e.g. ( Carter et al.|2011 Morales et al. 2009; 
|Irwin et al.|20TT| >). 
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